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ABSTRACT 

Some short-period exoplanets ( "hot Jupiters" ) are observed by their transits to have 
anomalously large radii. It has been suggested that these planets are in a resonance 
involving persistent misalignment and synchronous precession of their spin and orbital 
angular momenta — a Cassini state — and that the attendant tidal heating inflates the 
planet. We argue against this. Using explicit tidal integrations, we show that although 
an oblique Cassini state can dissipate many times the planet's rotational energy, the rate 
of dissipation must be much less than hypothesized. Dissipation causes the planetary 
spin to lie at an angle to the plane containing the orbital and total angular momenta. 
If dissipation is too rapid, this angle becomes so large that Cassini equilibrium is lost. 
A separate consideration limits the total energy that can be extracted from the orbit. 
The source of the torque on the orbit, either an oblique parent star or an inclined third 
body, aligns with the orbit as it absorbs the angular momentum shed by the planet. 
Alignment removes the orbital precession required by the Cassini state. In combination 
with observational bounds on the mass and semimajor axis of a possible second planet 
and with bounds on the stellar rotation and obliquity, these constraints make it very 
unlikely that obliquity tides can be the explanation for inflated hot Jupiters, especially 
HD 209458b. 

Subject headings: celestial mechanics, planetary systems, stars: individual (HD 209458), 
methods: analytical, methods: numerical 



1. Introduction 



Transit observations have exposed the limits of our understanding of the structure of short 
period extraso lar giant planets, or hot Jupiters. T he models of transiting planet HD 209458b 



constructed by iBodenheimer. Lin, fc Mardlingl (|200ll ) gave the first indication of a discrepancy be- 
tween the theoretically predicted and observed radii of the largest hot Jupiters. While standard 
evolutionary models for irradiated giant planets faithfully represent the majority of the 14 known 
transit cases, the conclusion persists that the models underpredict the radii of several hot Jupiters 
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by 10 - 20% (IBurrows et all l2006al : lArras fc Bildstenl 12009 : IFortnev et al.1 l2006al: iGuillotl 120050. 
At the large end of the size distribut ion, planets HP 2 09458b (|Knutson et al.ll2007l ). WASP-lb 
(jShporer et al,ll2006l ). and HAT-P-lb (jBakos et al.ll2006l ) present the most serious challenge to ex- 
pected rates of contraction. Resolution of the apparent radius problem will require simultaneous 
treatment of several outstanding issues, with competing effects. First, the problem is ameliorated 
somewhat by distinguishing the photospheric radius from the "transit radius," which is the impact 
parameter of photons that travel a chord of unit optical depth through the planetar y atmosphere. 
The l atter radius is larger by as much as ~ 5% for planets with low surface gravity (jBaraffe et al 



20031 ; IBurrows et al .1120031 . l2006al ). Secondly, significant enlargement ca n be achieved in two ways; 
either by tapping an additional heat source ( Bodenheimer et al] hoOll ) or by insulating the inte- 
rior with a supersolar abundance of metals in the atmosphere, thereby retarding its contraction 
(jBurrows et al.ll2006al ). Conversely, more rapid contraction may result if global redistribution of 
the i nsolation is inefficient, allo wing an increased cooling flux to escape via the planet's cold night 
side (jGuillot fc Showmanll2002l ). Finally, planet radii are greatly reduced by the inclusion of heavy 
elements in the interior, whether di spersed or as a solid core , for which a mass 10 — 40 times the 



mass of the Earth may be typical (jBodenheimer et al.l 120031 ; IFortnev et al.1 l2006al : IBurrows et al 
2006ah . 



Burrows et al.l (|2006al ) argue that, to within the observational uncertainties in radius and age, 



with no additional heating mechanism, the radii of all transiting hot Jupiters are consistent with 
evolutionary models if supersolar atmospheric metallicities are allowed. The authors rely upon 
heterogeneity in both atmospheric abundances and core mass. In particular, Burrows et al. fit 
coreless models with 10 times solar atmospheric abundances to WASP-lb (within the 1 a error 
bar), HAT-P-lb (within 1 a), and HP 209458b (within 2 a). A concern is that the inflated radii 
of these models require that the metals be sequestered in the atmosphere; supersolar abundance s 
throughout the planet might be expected to result in a net contraction (IFortnev et al.l l2006bl ). 
Observations of thermal emission from hot Jupiters also bear on the relevance of such models. 
Broadband detection s of infrar ed emission at secondary e clipse have been rep orted for HP 209458 b 



(jPeming et al.ll2005h . TrES-1 (jCharbonneau et al.ll2005T). and HP 1 89733b foeming et al 



Corresponding theor etical spectra were construc t ed by 



Fortnev et al.l (|2005l ). and IBurrows et al.l (|2005l . l2006bl ). All obtained fits to the emission mea 



Seager et al.1 (j2005h . Barman et al.l (|2005h 



2006) 



surements whi l e assum ing solar abundances. N one predicted metalli city as large as required by 
Burrows et al.l (|2006al ) for HP 209458b, although IFortnev et al.l (|2005l ) found marginally improved 
fits for TrES-1 using 3 — 5 times solar abundances. Interpretation should nonetheless be cautious, 
as the continua of the model spectra depend more strongly on the degree of horizontal heat trans- 
port than on metallicity. In fact there is no consensus among these works as to whether the data 
indicate that the reradiation of absorbed stellar flux occurs locally, is globally uniform, or is in- 
termediate between these limits. This points to another critical issue; insofar as the photospheric 
temperature is an indicator of the temperature near the base of the radiative zone, the observed 
day-night temperature contrast can constrain the redistribution of heat at depth. The contraction 
rates computed by Burrows et al. hinge on thorough redistribution and isothermal conditions along 
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the radiative-convective boundary, where the cooling flux is determined (jArras fc Bildstenl 120061 ) . 
Such an isoth erm bounding th e convective interior was seen in the numerical simulations of global 
circulation bv llro et al.l (|2005l ). However, if this deep isotherm does not exist, the n additional heat 



i ng is indeed required to offset the steeper temperature gradient on the night side. lHarrington et al 



(|2006l ) reported measurements of the phase curve of (non-transiting) hot Jupiter v Andromedae b 
that imply a day-night contrast perhaps exceeding 10 3 K. However, from simulations with a deep 
isotherm, Iro et al. predicted contrasts not larger than 600 K for HD 209458b (despite the fact that 
HD 209458b receives greater stellar flux than v And b). As a final point, the first spectroscopic 
measurements o f thermal emission h ave come with thei r own surprises. The mid-infrared spectra 



of HD 189733b (iGrillmair et al.ll2007h and HD 209458b dRichardson et al. 



2007 



Swain et al.ll20071 ) 



are consistent with a Rayleigh-Jeans continuum, lacking the anticipated absorption bands blueward 
of 9 /Um associated with water and methane. These observations of hot Jupiter atmospheres will 
need to be understood theoretically before a conclusion can be reached regarding heat retention. 
A continual heating mechanism may yet play a principal role. 

Augmentation of a Jovian planet's intrinsic luminosity causes an adjustment to a new hydro- 
static equilibrium with an expanded radius. The estimated heating rates necessary to account for 
the inflated hot Jupiters are in the range 10 26 — 10 28 ergs _1 , which is < 5% of the intercepted 
stellar luminosity and is a few orders of magni t ude larger than the intri nsic cooling luminosity 



(Burrows et al 



2006a 



Bodenheimer et al 



2003 



2001 



Baraffe et al 



2003] ) . Th is energy may be 



drawn from either the stellar insolation or the orbit. iGuillot fe Showman! (|2002l ) pointed out that 
one-sided irradiation should drive an atmospheric heat engine that converts the insolation into bulk 
flow kinetic energy. Inflation would require that a fraction of the kinetic energy be dissipated at 
depths where the radiative timescale is sufficiently long. Numerical global circulation models that 
address this mechanism have yet to conv erge on several key factors, including characteristic wind 



Showman Sz Guillotj ( 2002 ) found substantial downward 



Burkert et al 



(j2005h bund little vertical tra n sport and 



velocities and global flow patterns. While 
transport of kinetic energy in simulation, 
near complete dissipation at shallow depths. A further puzzle is that lHarrington et al.l (|2006l ) ob- 
served the phase curve of v And b to be consistent with an 11° lag between full ph a se an d the 
substellar point. A super rotational wind like that predicted by IShowman &: Guillotl (|2002l ) and 
Cooper &; Showman! (|2005l ) for HD 209458b, with a 60° advance, is strongly disfavored by the data. 
Presuming that dissipation does occur at sufficient depth to affect the radius, one has to explain 
why this mechanism — or any other that appeals to characterstics common to all hot Jupiters — does 
not inflate all currently observed transiting planets. A possible answer involves competition among 
conflicting factors: for example, auxiliary heating offset by variable core masses. 

An alternative is that tidal forces may heat the planet at the expense of the orbit. Tidal 
dissipation naturally depends upon variable orbital and perhaps other parameters. The action 
of tides is to align the spin and orbital angular momenta, synchronize the rotational and orbital 
periods, and circularize the orbit. The efficiency with which hot Jupiters are thought to dissipate 
tidal oscillations implies that alignment and synchronization occur on 10 5 — 10 6 yr timescales and 
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circularization occurs on 10 8 — 10 9 yr timescales. The eccentricity and/or obliquity must be excited 
continually in order for tides to persist for several Gyr. Unlike the case of insolation-driven heating, 
ongoing tidal heati n g req uires special orbital conf i gurati ons, potentially involving a third body. 
Bodenheimer et al.l (|200ll ) and iBodenheimer et al.l (|2003l ) have suggested that eccentricity could 
be forced by a companion planet, and they have estimated that a sustained eccentricity e > 0.03 
would be sufficient to explain the radius of HD 209458b . Recent radial vel ocity measurements 
(jLaughlin et al.ll2005al ) and timings of secondary eclipse (jDeming et al.ll2005l ) have tightened the 
constraints on the eccentricity of this system, but because of degeneracy with the longitude of 
pe riastron, cannot y e t rule out e = 0.03 with 95% confidence. An innovative possibility, put forth 
by IWinn &: Holmanl (|2005l ). is that a planet could occupy a so-called Cassini state in which its 
obliquity, rather than its eccentricity, is forced. 

G. D. Cassini published three empirical laws in 1693 describing the spin state of the moon. 
Cassini's laws state that (1) the moon's rotational and orbital periods are synchronous, (2) its 
equatorial plane maintains constant inclination to the plane of the ecliptic, and (3) its spin axis 
remains coplanar with the normal to its orbit and the normal to the ecliptic. The first law is an 
outcome of tidal evolution for bodies with a permanent quadrupolar distortion, like the moon, or 
for fluid bodies in circular, stationary orbits ( Murray Sz Dermottlll999l ). The second and third laws 
characterize a Cassini state, which is an equilibrium for the motion of the spin axis of a body in 
a uniformly precessing Keplerian orbit. Up to four Cassini states may exist for a given system, 
each with a preferred obliquity. The moon occupies Cassini state 2 (see numbering scheme in £j2|) 
with obliquity 9 ~ 6? 7. Resonance between the precession frequency of the orbit and the natural 
precession frequency of the spin determines the values of the Cassini obliquities. The obliquity of a 
dissipative body is driven toward a Cassini state, which is a quasi-equilibrium in the sense that the 
body continues to experience an obliquity tide so long as 6 is nonzero. In this paper we determine 
the effect of the attendant obliquity tide on the evolution and survival of Cassini states for hot 
Jupiters. 



Levrard et al.1 (|2007l ) have recently argued that the probability of a planet becoming trapped 



in an oblique Cassini state is rather small. Such probabilities depend upon the early evolution of 
the star and perhaps of its protostellar disk. Independently of initial conditions, however, there are 
reasons to doubt whether a hot Jupiter can remain in Cassini st ate 2 while experienc ing ongoing 
tidal heating sufficient to affect the planet's radius measurably. iLevrard et al.l (120071 ) themselves 
have pointed out that maintenance of the oblique Cassini states requires a balance between the 
torque associated with dissipation in the planet and a nondissipative torque on the orbit due either 
to an oblate and oblique stellar primary or to a third body. Since the nondissipative torque on the 
orbit is likely to be small, we argue that the Cassini state can be maintained only by assuming a 



dissip ation rate in the planet that is too small to inflate its radius significantly. Also, ILevrard et al 



(|2007l ) considered only the balance of torques tending to alter the obliquity. Using a specific tidal 
model, we show that a stronger constraint results from requiring the spin to equilibrate in the 
azimuthal direction. 
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Our investigation, however, was sti mulated not by th e poin ts above but by concerns about the 
orbital angular momentum. As noted by I Winn &; Holmanl (120051) . the suggested tidal i nput of > 4x 
10 26 erg s _1 required to explain HD 209458b's present radius ( Bodenheimer et al.ll2003l ). if continued 
for 5Gyr, is twice the present orbital energy of the planet. Hence, even if a somewhat smaller 
tidal power were acceptable or the planet suffered less dissipation in the past, the orbital energy 
must have chang e d by of order its present value, with a large loss of orbital angular momentum. 
Winn fe Holmanl ()2005l ) have not explained where the angular momentum might have gone. It is 
not likely to have been absorbed by the star, for two reasons. First, the moment of inertia of the 
star is somewhat smaller than that of the orbit (£ )3.3p . so that, notwithstanding magnetic braking 
by a stellar wind, stars hosting inspiraling hot Jupiters might be expected to rotate more rapidly 
than other stars of similar age and spectral type, contrary to the evidence of Figure [3] in this paper. 
Second, the timescale for dissipating the planetary tide in the star probably exceeds its age ( ^3.3|) . 
which limits the amount of angular momentum that the star can absorb. Without dissipation, the 
star can change only the direction but not the magnitude of its spin angular momentum^ Thus, the 
amount of angular momentum that the star can absorb from the orbit is limited by the star's initial 
obliquity as well as its moment of inertia. In four cases where the spectroscopic transit (Rossiter- 
McLaughlin effect) has been well measured, the implied angle (A) between the p rojections o f the 
stellar spin and the orbital angular momentum onto the plane of th e sky is small rtWinnll2006l . and 



references therein); in particular, A = -4? 4 ± 1?4 for HD 209458 (jwinn et alJbood ). Unless the 



stellar s pin has a large comp onent along the line of sight — which has been ruled out for HD 189733 
at least (jWinn et al.ll2006bl ) — the smallness of A indicates that if maintenance of the Cassini state 
relies upon angular momentum changes in the star, then the future life expectancy of that state is 
much less than the present age of the system. 

The plan of this paper is as follows. $2] presents a pedagogical review of Cassini states in 
general nondissipative, Hamiltonian systems. discusses causes of the precession of the orbit 
in more specificity: either an oblate and oblique primary, or a third body. Tidal dissipation is 
introduced, and it is argued from empirical constraints on tidal quality factors that dissipation 
in the planet is more important than dissipation in the star for this problem. SJH presents and 
interprets the governing equa tions used to describe the tida l evolu tion of the planetary spin and 
orbit, which are adopted from lEggleton &; Kiseleva-Eggletonl (|200ll ) with small notational changes. 
In fJH we present exemplary numerical results and interpret them using analytic approximations. 
We explore how Cassini state 2 "breaks" due to the imbalance of dissipative and non-dissipative 



1 In the equilibrium tide, the tidal bulge raised in the star by the planet points directly at the planet, since the 
phase lag between the bulge and the perturber is proportional to the dissipation rate. Hence the planet exerts 
no torque on the tidal bulge, while any torque exerted on the star's rotational oblateness is perpendicular to the 
spin axis and therefore changes only the direction of the spin. Alternatively — and the following argument applies 
to dynamical as well as equilibrium tides — we may apply Kelvin's Circulation Theorem to a closed contour drawn 
along the rotational equator. The circulation 2nR±v TO t is conserved in a nondissipative star, since (V X v) -V S is 
conserved following the fluid and we may take the curve to be the boundary of a surface of constant specific entropy 
S. A contradiction would therefore result if the equatorial velocity were to evolve secularly without dissipation. 
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torques, and we find a strong constraint on the planetary dissipation rate that renders obliquity 
tides unimportant as an internal heat source. In £j6l we summarize our main results and use them 
to quantify the objections raised above to Cassini states as an explanation for the anomalously 
large radius of HD 209458b in particular. 



Cassini states 



The spin dynamics of a body in a precessing Keplerian orbit are analogous to orbital dynamics 
in a rotating 2D potential. In both systems a corotating frame is found in which the stationary 
points of the effective potential are equilibria. The spin counterparts to Lagrange points are Cassini 
states. In this section we review the theory of Cassini states. 



The modern study of spin equilibria w as begun by I Colombo! (119661 ). who generalized Cassini's 
laws for axisymmetric bodies. iPeald (|1969l ) first enumerated the Cassini states and generalized to 
the non-axisymmetric case. Initial applicat ions of the theory were made primarily to Mercury and 



the m oon, but also to Iapetus and Triton (lColombdll966l : iPeale 



1969 



1974 



Beletskiilll972l : I Ward 



19751 ). Evolutio n of the ob li quity of M ars due to geophysical c hange s in the planet's oblateness 
was studied by I Ward et al.l (|1979l ) and iHenrard &: Murigandd (j 19871 ) . The latter work contains 
a particularly illuminati ng developme nt of the Cassini form alism. Recent applications to Saturn 
dWard fc Hamiltonll2004h and Jupiter dWard fc Canirokood ) suggest that the primordial obliquities 
of many solar system bodies may have been alt ered by trapping i n Cassini resonances . Our present 



concern is with application to hot Jupiters (jWinn fc Holmanl 120051 ; iLevrard et al.l 120071 ) . The 
essential features of Cassini states may be deriv ed by s pecia lizing to the case of a dissipationless, 
axisymmetric body. For a general treatment see IPeald ()1974l ). We discuss the critical role of tidal 
dissipation in §4.21 



Consider an axisymmetric planet with principal moments of inertia C > A = B in orbit around 
a star of mass M*. Let the total angular momentum of the system (including other bodies) be 
denoted J, pointing along the fixed unit vector J. This axis and the invariable plane perpendicular 
to J are at rest in an inertial frame. Let the orbital plane be inclined to the invariable plane by 
angle /, and allow the node to precess uniformly around J with angular frequency g due to an 
unspecified constant torque. In the two-body problem the nodal precession is dominated by the 
stellar quadrupole, but torques from additional bodies or the circumstellar disk may dominate when 
they are present. The orbit normal h remains fixed in a frame rotating with angular velocity g = gJ 
with respect to the inertial frame. In cases of interest the node regresses (g < 0). Define a Cartesian 
coordinate system XYZ attached to the rotating frame as in Figure [IJ with the Z-axis along the 
orbit normal (Z = h) and the X-axis in the direction of the ascending node. The orientation of 
the planet relative to the rotating frame is specified by the usual Euler angles (6,<p,tp), which give 
the positions of the planet's axis of symmetry z and the arbitrary orthogonal axes x and y frozen 
into the body. The equatorial plane intersects the orbital plane at an obliquity 6 = cos _1 (z • h). 
This node precesses around h with angular frequency <j). The planet rotates with respect to the 
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inertial frame with instantaneous angular velocity u, which is the sum of the spin angular velocity 
= ijjz and the angular velocity associated with the precession of the symmetry axis. If the 
z-axis is at rest in the rotating frame then cj = CI + g. The planet's total spin angular momentum 
S = Cu z z+Bcu y y+ Au x x can be approximated by principal axis rotation if the precessional angular 
frequency is small compared to f2, in which case only the z component of S is retained. The further 
approximation lo z ~ O is often made for simplification, and we do so here. In Figure [T] we have 
implicitly assumed principal axis rotation by taking S = Cfi. The Hamiltonian Wrf governing the 
spin of the planet in the rotating frame is related to the Hamiltonian TCif in the inertial frame by 
T~Lrf = T^if — S • g. This is simply the relation between Jacobi's integral of motion in a rotating 
potential and the total energy. TCip is the total energy associated with the planet's orientation, and 
Wrf is an integral of motion for constant g and in the absence of dissipation. Hereafter we drop 
the subscript RF. The Hamiltonian in the rotating frame 



c2 



S-g 



(1) 



includes the orientation-dependent component V of the planet's gravitational potential. This is 
contributed by spin-orbit coupling, while spin-spin coupling by the stellar quadrupole is negligible. 
Expanded in spherical harmonics and truncated at second order, 



V 



GM*(C-A) 



;(r-S) 



(2) 



where f is the instantaneous direction of the planet with respect to the star and r its distance. V 
has a short-term variation over the orbit and a long-term variation with the spin precession. We 
may average V over an orbital period while holding S constant to obtain the mean potential on 
timescales much greater than the orbital period P. For an orbit with mean motion n = 2tt/P and 
eccentricity e, 



1 

P 



Vdt 



n 2 (C-A) 



(S • hf 



(3) 



, 2(1 -e 2 ) 3 / 2 

which depends only on S since h is fixed in the rotating frame. The system described by (pQ) 
and ([3]) has 2 degrees of freedom. Its state may be specified by the generalized coordinates (ip, (ft) 
and conjugate momenta (S, S cos 6). It is immediately apparent that S is conserved due to the 
axisymmetry of the planet. We define the constants 

2C 



a = 



2Cn(l - e 2 ) 3 / 2 



S' 



a 
6' 



The former characterizes the natural precession frequency of the spin. The latter is the part of the 
Hamiltonian composed of variable terms: 



H' = -g(S-J) 



f(S-h) 2 . 



(4) 



Viewed in terms of the projection S = (X, Y, Z) onto the r otating axes, e quation @) describes a 
family of parabolic cylinders opening in the +Y direction ([Colombo! I1966I ). Solution trajectories 
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for S are curves of intersection between the parabolic cylinders and the unit sphere. The canonical 
equations of motion cp = &H/d{S cos 9) and d(S cos 6)/ dt = —dTt/dcp become 

— = g(S X J) + a(S • h)(S X h). (5) 

The equilibrium solutions of ([5]) are Cassini states. Constant obliquity requires that the spin axis 
be coplanar with J and the orbit normal (i.e. (j) = or tt), because otherwise S X J and S X h 
point in different directions. This statement is Cassini's third law. Constant (f> requires that the 
equilibrium obliquity 9 satisfy a relation 

= g sin(# — I) + a cos 9 sin 9 (6) 

that has either two or four roots for given values of g/a and /. By convention 9 is negative when 
<j) = and positive when <j) = ir. Equation ([6]) is identical to that obtained by finding the stationary 
points of 7i ! on the unit sphere, which are points of tangency between the unit sphere and parabolic 
cylinders. 

Figure [2] illustrates the equilibria admitted by equation ([6]) for / = 0.1 (left panel) and / = 1.0 
(right panel). In the limit that the orbit precesses slowly (|<7/a| <C 1) the term derived from V 
dominates the Hamiltonian. The mean gravitational potential is minimized when the spin is aligned 
or antialigned with the orbit normal (states 1 and 3, respectively), and maximized when the spin 
tips into the plane of the orbit, pointing in the +Y direction (state 2) or the —Y direction (state 
4). States 1, 2, and 3 are stable against small displacements, which result in libration around the 
equilibrium point. The torque S X g is responsible for libration within a domain around state 2, 
which would otherwise not be stable. State 4 is a saddle point of TC', lying on the trajectory that 
separates the domains of the other three states (i.e. the separatrix), and is therefore unstable. 
As \g/a\ increases the domain of state 2 expands while that of state 1 contracts, until a critical 
point is reached at which states 1 and 4 merge and vanish. In the limit of a rapidly precessing 
orbit (|g/a| S> 1) the planet responds to the time-average position of the orbit normal, namely J. 
H! is then maximized for spin aligned with J (state 2) and minimized for spin antialigned with J 
(state 3). Comparison of the left and right panels illustrates that, at larger inclinations, the domain 
around state 2 is broader because it lies closer to J. 



Winn Sz Holmanl (120051 ) pointed out that Cassini state 2 is the most favorable configuration 
wherein a planet could maintain a significant obliquity. The obliquity of state 1 is large only for 
a coincidental similarity of g and a (near the critical point), which are unrelated parameters. As 
for state 3, we show in §4.11 that no such equilibrium exists in the presence of a dissipative tidal 
torque. 



3. Tidal interactions 



As we have just seen, a torque must be exerted on the planetary orbit to keep it in Cassini 
state 2. This torque could be tidal, involving the oblateness and obliquity of the central star, or 
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it could be due to a third body in an inclined orbit. We examine both possibilities below. Tidal 
dissipation must occur in the planet at some level if it remains in Cassini state 2. We compare 
the relative importance of dissipation in the planet and in the star for the evolution of the orbital 
energy and angular momentum. 



3.1. Orbit precessed by stellar oblateness 



Given dissipation in a body, the spin will move towards alignment with the orbit normal. 
However, given reasonable values for the star's energy dissipation rate, the tide of a hot Jupiter is 
too weak to cause the angle between the star's spin and its orbit to change substantially (the time 
scale for such change is s imilar to equation (1131) below). Therefore, measuremen ts of alignment 



among extrasolar planets (jQueloz et al 



2000 



Winn et al.ll2005l : IWinn et al.ll2006al ) are understood 



to be faithful indicators of the primordial inclination. A system that is not perfectly aligned 
will undergo precession analogous to that of the Earth's spin as the Earth's rotational bulge is 
torqued by both the moon and the sun. As an extrasolar planet causes its host star to precess, 
the planet's orbit also precesses, conserving total angular momentum. By this mechanism, the 
orbits of hot Jupiters m ay precess at rates that are fast enough to be observable by transit timing 
(jMiralda-Escudel 12002 ). In such systems, the star's spin angular momentum S* = C+fl* has a 
magnitude comparable to, but somewhat smaller than, that of the orbital angular momentum L, 
and the planet's spin angular momentum S p = C p ft p has negligible magnitude. The total angular 
momentum J = L + S* + S p is located between the stellar spin and orbit normal, but it lies slightly 
outside of the plane defined by L and S* unless S p lies in that plane. If the time average of S p is 
kept out of that plane, then the misalignment between stellar spin and the orbit normal can evolve, 
even without dissipation in the star (see $5|). 



3.2. Orbit precessed by third body 

Another source of orbital precession is a third body. There are two known transiting planets 
in binaries (HD 189733b & HAT-P-lb), but the stellar companions would cause planetary orbital 
precession on a timescale much longer than the host star's stellar oblateness would, so their effect 
is negligible. No system with a transiting planet has a known second planet. 

Suppose there is an undetected second planet on an exterior orbit. It will perturb the eccen- 
tricity of the transiting planet on a short timescale, which will be continually damped by solar 
tides, which would also heat the planet. However, as the transiting planet's orbit shrinks due to 
this continual excitation and damping of eccentricity, the rate \g\ at which it is precessed by the 
exterior planet will decrease. Also, as the orbit shrinks, a increases. If the planet started with 
\g\/a > (|g|/a)criti Cassini state 2 was the only equilibrium available. As \g\/a slowly drops, the 
spin axis will tip over as \g\/a crosses (\g\/a) cr it and the spin remains in state 2 (see Fig. [2]). There- 
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fore a mixture of eccentricity and obliquity heating is possible, if Cassini states are stable despite 
dissipation. 

The amount of angular momentum that the second planet can absorb from the first as its orbit 
shrinks can constrain the parameters of the second planet. Let L and L 2 be the orbital angular 
momentum of the transiting planet and the second planet, respectively. Suppose that L and L 2 
are the dominant angular momenta of the problem such that, by the law of cosines, 



J 2 = L z + Li.+ 2LL 2 cos i re i, (7) 

where i re i is the relative inclination of the two planets' orbits. The semimajor axis and eccentricity 
of the second planet are both conserved on secular timescales if the semimajor axis ratio is large, 
which can be shown by taking th e potential of the second planet to quadrupole order and integrating 



over the orbits of both planets (IFord et al.l 12000). Therefore, the only way the orbit of a second 



planet can absorb angular momentum from the dissipating planet is by reorienting. As the planets 
approach coplanarity, the orbital precession period becomes too long for the dissipating planet to 
remain in Cassini state 2. The obliquity damps to zero and loss of energy and angular momentum 
halts. For the final condition of i re \ = 0, we therefore have J = Lf + L 2 , where Lf is the final 
angular momentum of the dissipating planet. Constancy of J gives: 

L f = {Lf + L 2 2 + 2L i L 2 cosi rchi ) 1 / 2 - L 2 , (8) 

where Li and i Te \ j are the initial angular momentum of the dissipating planet and the initial relative 
inclination, respectively. Solving for the angular momentum of the second planet gives 

1 Lf - L 2 f 

L 2 = f - . (9) 

2 Lf - Li cosi re i,i 

For any solution to exist, we must have i r el,i > &'rei,crit = cos~ 1 (L/-/Lj). This critical value corre- 
sponds to the case in which L 2 S> L and essentially all of the reorienting is done by L. 

If we suppose that a particular planet is currently in Cassini state 2, with its semimajor axis 
still shrinking, then L QOW > Lf. If a second planet is the cause of orbital precession, its angular 
momentum must be greater than L 2 of equation (J9]), taking Lf — > L now . This constraint may be 
very restrictive, for a presumed value of i re i,i, if we can estimate Lj. Tidal dissipation is a strong 
function of semi-major axis; from the equations of $1] one may derive an approximate relation 
L oc — L~ 12 . If the rate of angular momentum evolution L is observed or is theoretically inferred 
(for instance, by an inflated radius), this relation can be integrated back in time for the age of the 
system to obtain Li. Auxiliary theoretical arguments may constrain i re i,i, so the angular momentum 
L 2 of an undetected second planet may be given a lower bound, if it is forcing obliquity tides. The 
example of HD 209458b will be examined in £JBJ 
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3.3. Dissipation 



The nature of tidal dissipation in Jovian planets and late-type main-sequence stars remains 
a subject of research. Such bodies are largely fluid, and the tidal Reynolds number based on the 
true microscopic viscosity is en ormous. Th erefore, most theories for the dissipation invoke resonant 
coup l ing to ste llar convection (IZahnlll966l ). or to modes of oscillation (IZahnlll970l : lOgilvie &: Lin 
20041 : IWull2005l ). Where these mechanisms can be worked out with confidence, they often fall short 
of explaining observationally inferred rates of circularization or synchronization. Other mechanisms 
show promise but require difficult and uncertain calcul ations. Therefore, in the p resent work, we rely 
upon empirical scalings from observations. Following iGoldreich &: Soterl (119661 . henceforth GS66), 
we parametrize the dissipation by a tidal quality factor Q defined as the ratio of the maximum 
stored energy in the tidal distortion of the star to the energy dissipated per radian of the tidal 
cycle. The "equilibrium tide" approximation is used, as if the tide were static, so the stored energy 
is potential. In this formulation, the phase lag between the applied tidal potential and the response 
is (2Q) _1 . The torque on a nonrotating body of mass M\ and radius R\ due to a companion of 
mass M2 at separation r is 

_ 3ki GM$R\ 



Qi 



(10) 



The quantity k\ is the apsidal- motion constant (which is half the Love number), which GS66 effec- 
tively replace by 3/4 as appropriate for a body of uniform density — a reasonable approximation for 
a solid planet or moon. The structures of solar-type stars and Jovian planets are better approxi- 
mated by polytropes of indices n* = 3 and n p = 1, respectively, for which (A;*, k p ) ~ (0.014, 0.26). 
In all tidal dissipation rates, however, Q and k appear in the same combination that one sees in 
equation ()10p above, so it is convenient to frame the discussion in terms of an effective quality 
factor! Q' = Qx 3 /(4k). 

Values for the quality factor can be inferred from orbital properties of appropriate binary stars, 
extrasolar planets, and satellites of Solar System planets. Since it is simply a parametrization of 
ill-understood dissipative processes, Q' may depen d upon tidal period, and also on subtle details 
of internal structure. As an example of the latter, lOgilvie fc Linl ([2004J) state that their inertial- 
wave mechanism depends upon the existence of a rocky planetary core; there is evidence from the 
measured radii of transiting planets that their cores vary considerably in size and mass. 

With these caveats in mind, we review some of the observational inferences regarding Q'. GS66 
estimated that the mean-motion resonances of Io, Europa, and Ganymede require Jupiter's quality 
factor to lie in the range 10 5 < Q'j < 10 6 , unless that resonance is accidental or short-lived. The 
relevant tidal period in this case is 6.5 hr, one half the synodic period between Jupiter's rotation 
and Io's revolution. For spectroscopic binary stars, tidal circularization is easier to recognize 
than tidal synchronization, since orbital eccentricity is measurable from the radial-velocity curve. 



2 In GS66's original notation, Q and Q' differ by a term involving elasticity, which can be significant for solid 
planets but is presumed to be negligible in our case. Our Q' reduces to their Q for a uniform-density body. 
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Mathieu et al.l (|1992l ) and lMathieu et al.l (|2004l . and references therein) find that the circularization 
period (-P c irc) separating circular from measurably eccentric binary orbits in stellar clusters increases 
with cluster age beyond 1 Gyr, indicating that significant tidal dissipation occurs on the main 
sequence. For M67 and NGC 188, for example, P c i rc = 12.4 d and P c i rc = 15 d are quoted, with 
corresponding cluster ages of 5 and 7 Gyr, respectively, whereas clusters younger than 1 Gyr are said 
to be consistent with P c i rc ~ 8 d, perhaps reflecting rapid circularization during protostellar phases 



(jZahn fc Bouchetlll989l ). To relate such observations to tidal dissipation rates, it is conventional 
to suppose that P c irc is the period at which eccentricity has decayed by three e-foldings. Adapting 
GS66's formula for the damping rate due to star 1 alone, we have 



dine 
dt 



171 M 2 
16Q' Mi 



a j 



n. 



(11) 



in which a is the binary semimajor axis, n = (GMtot/a 3 ) 1 ^ 2 is the mean motion of the orbit, and 
synchronous rotation has been assumed. Applying this to a binary consisting of two sunlike stars 
and summing the dissipation rates in both stars, we estimate from M67 that Q'q ~ 3.5 x 10 5 , and 
from NGC 188 that Q' Q » 2.2 x 10 5 . 

The fact that the orbits of most extrasolar planets with periods less than ~ 5 days and masses 
> Mj are essentially circular provides another constraint on Q' if we assume that this is due to 
tidal dissipation in the planet o Again taking \dlne/dt\~ l = t/3 at P c i rc and using t = 5 Gyr as a 
representative age, one has 



Q'«3x 10 5 



5d 



-13/3 



-2/3 



(Mj_ 



Alternatively, circularization might be due to dissipation in the star. However, equation (jllj) 
predicts that dlne/dt is dominated by whichever body has the smaller value of Q'M 2 /B? . This is 
probably the planet, as discussed below. 

It is reassuring that these various estimates for Q' in Jupiter, extrasolar planets, and sunlike 
stars are comparable; in fact, one might boldly conclude from the evidence above that there is 
a universal value Q' ~ 3 x 10 5 . Many theories of tidal dissipation, however, predict that Q' 
should scale with tidal period. For equilibrium tides subject to a frequency-independent effective 
viscosity (perhaps associated with turbulent convection), it can be shown that Q' oc P. With 
this scaling, the values estimated for Jupiter from its Galil ean moons should be scaled u p by an 
order of magnitude for extrasolar planets in ~ 4 d orbits. iGoldreich fc Nicholson! (119771 ) argued 
that when the tidal period is less than the turnover time of the energy-bearing convective eddies, 
the effective viscosity they provide should be suppressed by a factor ~ (Ptide/Peddy) 2 ; this would 



3 At the time of writing, based on the catalog at http://vo.obspm.fr/exoplanetes/encyclo , we find P c irc > 5d 



GJ 436b is the only planet with P < 5d that has an eccentricity securely above 0.1, which iManess et al 
attribute to either ongoing eccentricity pumping or an unusually high Q for this rather low mass (M sin i « 
planet. 



( 20071 ) 

23M ffi ) 
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lead to a quality factor decreasing with increasing period as Q' oc P _1 . iGoldman Mazehl ()199lh 
suggested that this scaling fits the observed trend of binary circularization periods with age better 
than an increasing or constant Q' . However, if the scaling Q' oc P _1 results from a suppression 
of convective dissipation at short per iods, then it is very un li kely t hat convection can account for 
the inferred Q' of Jupiter a t all, as iGoldreich &: Nicholson! ([19770 point out. The inertial-wave 



mechanism of lOgilvie fe Linl (|2004l ). on the other hand, apparently predicts a quality factor that is 
independent of tidal frequency for a fixed ratio between the tidal and rotation frequencies. 

In the case of hot Jupiters, dissipation in the planet is likely to be more important than 
dissipation in the star. On the one hand, it follows from equation (|10p that the ratio of the 
synchronization torques on the spins of the two bodies is 



r, 



91 



n 2 R 



fiR* 



(12) 



where pi = 3Mj/4-7rP 3 is the average density of body i. The same ratio determines the relative 
importance of the two bodies for eccentricity damping. The densities of the star and planet are 
typically comparable, and the mass ratio M p : ~ 10~ 3 , so the second factor on the right hand 
side above is < (Mp/M*) 1 / 3 ~ 0.1. So the torque on the planet is greater unless the star has a 
substantially smaller Q' . On the other hand, the synchronization timescale due to the torque on 
the star is approximately 



t 



(M p a 2 + C* 



sync,* 



nh — fij 



15 



3 x 10 5 



P 

Id 



13/3 



M 5/3 

Pq) 



Mj 



Gyr, (13) 



in which C* ~ 0.06M+R 2 is the stellar moment of inertia. Note that the stellar and orbital 
angular velocities are differenced as vectors, to allow for possible stellar obliquity. In the final 
numerical form, we have neglected f2* compared to nh and taken C+ ~ M p a 2 : in fact C+/(M p a 2 ) ~ 
0.6(P/4d)~ 4 / 3 for a sunlike star and a Jupiter-mass planet. Thus, the stellar spin will not have 
been substantially altered by its tidal interaction with the planet unless Q'+ < 10 5 , contrary to 
the evidence of binary stars. Of course, Q'+ may depend upon the companion. Theory suggests 
that Q* may be gre ater under the tidal influence of a planetary rather than a stellar companion 
(|Ogilvie Linll2007h . 



Moreover, the spectroscopic (v sin i)+ values of exoplanetary hosts indicate that the stars rotate 
subsynchronously. To emphasize this point, in Figure [3] we plot the values of (vs'mi)+, relative to 
the synchronous value, for the known transiting planets. The value of sin i should be very close to 1 
if ther e is approxim ate alignment between stellar spin and planet orbit, as has been observed in four 
cases (|Winnll2006l . and references therein). All of the stars rotate much slower than their planet's 
orbit. In fact, their spin rates are largely consistent with a standard Skumanich law — an empirical 
relation v oc t -0 5 that describes the braking with age t of late-type dwarf stars due to the stellar 



wind (ISkumanichlll9720 . A suitable reference for calibrating the Skumanich law is lValenti &: Fischer 
(|2005|)'s spectroscopic catalogue of 1040 nearby F, G, and K stars that have been targeted by planet 
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searches. We calibrate separately for six mass bins, the first with 0.7 < M*/Mq < 0.9, and the 
remaining five bins spaced between 0.9 — 1.4Mq, each with width 0.1 Mq. In all bins a better fit 
to the data can be obtained with a power law steeper than —0.5, but we ignore this systematic and 
depict standard Skumanich laws for simplicity. Scatter around the fits is « 10%. Only HD 209458 
appears to be rotating more rapidly than would be predicted by magnetic braking alone (and only 
marginally so, given the scatter and systematic uncertainties). The tidal torque these stars receive 
from their planets is apparently weak compared to the torque from their magnetospheres. Stars 
HAT-P-1 and TrES-1 are the relatively slow rotators of the group; however, note that on this 
figure, the Sun would also be rotating slowly by a factor of 2. 

We conclude that for the purposes of the present paper, tidal dissipation in the host star can 
be neglected. 



4. Spin and orbit equations of motion 



Here we introduce the equations of lEggleton fc Kiseleva-Eggletonl (120011 ) for the spin and orbit 
evolution of two fluid bodies. The equations of motion are for the orbital elements in the secular 
approximation. That is, the Newtonian equations of motion are integrated over one orbit, and 
what remains is time derivatives for the orbital elements. This formulation allows us to perform 
integrations quickly that represent dynamical evolution over billions of years. The parameters of 
energy dissipation in this model are based on the ass umption that the ener gy loss rate is proportional 
to the time rate of change of the quadrupole tensor (jEggleton et al.lll998l ). A prolate bulge is raised 
on each body by the tidal potential of its companion. If the rotation is not synchronized, energy 
dissipation drags the bulge away from the instantaneous direction of the companion. The familiar 
consequences are that the spins of the bodies synchronize with the orbit and that the orbit's 
eccentricity is damped as its energy is sapped at nearly constant angular momentum. Another 
effect quantified by this model is the spins coming into alignment with the orbit normal, which 
is the key piece of technology we need to investigate the relationship of tidal synchronization to 
Cassini states. 

We shall specialize to the limit of zero orbital eccentricity, for a few reasons. First, the 
eccentricity will damp to zero on a timescale short compared to the age of the system. Second, the 



effect of eccentricity damp ing in hot Jupiters has been pursued by other authors (IBodenheimer et al 



200ll : iLevrard et al.ll2007l ) , so this approximation more clearly focuses attention on the contribution 
of obliquity to tidal dissipation. Finall y, this approximation considerably simplifies the equations 
of lEggleton fc Kiseleva-Eggletonl (|200ll ); for instance, the only included effect of general relativity 
is to cause the apse to precess, which has no physical effect when e = 0. 

Let the two masses equal M* and M p , of total mass M, and of reduced mass fi = M+M p /M . 
Let h denote the orbital angular momentum per unit reduced mass, of magnitude h and direction 
h, and let e and q be reference directions in the plane of the orbit such that e X q = h. Let 
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the vector spins of the two bodies be fi* and f2 p with rotational moments of inertia C* and C p . 
Therefore, with our specialization to e = 0, there are 9 dynamical variables: the 9 components of 
the vectors h, f2*, and fi p . The equations conserve the total angular momentum. The equations 
of motion are: 



a 



ldh 

h~dt 



Q>* + y p )e - (** + Xp)^ - (W* + W p )h 
^(-^e + ^q + Wih), 



where 
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(14) 
(15) 

(16) 
(17) 
(18) 



tFi\ n 

and i is either p or *. i?j is the radius of each body, a = h?/(GM) is the semi-major axis, 
n = (GM) 2 h~ s is the mean motion, ki is the apsidal motion constant, and tpi is a tidal friction 
timescale: 



tFi 



4QJ 



Mf 



-n 



(19) 



9 \RiJ fiM 

The integrations in ^5]hold t he product Q' i n~ 1 cons tant during the evolution, as recommended by 
the equilibrium tide model of lEggleton et al.1 (|1998l ). 



4.1. Equilibration 

In equations (|16|) and (|17p . the first term on the right hand side is a non-dissipative torque 
that gives rise to nodal precession of the spins and orbit. The second term in these equations and 
the term on the right hand side of equation (|18p are due to tidal dissipation. Since dissipation in 
the star is negligible, we may set tp* — > oo, and specialize the discussion to dissipative torques on 
the planet. 

In equilibrium the magnitude £l p of the planetary spin angular velocity is constant. Therefore 

<r • n„) = o, (20) 

where T is the total torque on the planetary spin, and the angular brackets denote an average 
over the orbital mean motion. It suffices to replace T by the tidal (i.e. dissipative) component, 
since the precessional (i.e. nondissipative) part already satisfies (|20p . In the equilibrium tide, the 
instantaneous tidal torque takes the form 

T f = -Kfx (r X u), (21) 
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where u = nh — ft p and r is a unit vector from the star toward the planet. The logic behind 
equation (|2ip is that if the tidal bulge is a prolate distortion with major axis direction b, then 
Tt oc b X f . Absent dissipation, 6 = f so that the torque vanishes. Dissipation gives the bulge a 
phase lag of magnitude ~ (2Q)" 1 with a direction determined by the angular velocity of the tide seen 
in a frame spinning with the planet. This angular velocity is u, so that b ~ f+(2Q)~ 1 (f X u)/|f X u|. 
Comparison with the equations of motion shows that k = jj 1 h(ntp p )~ 1 . Substituting equation (|2ip 
into equation (f20l) gives 

ft p -(u-u-(ff)) = 0, (22) 

since the precession of fi p &: u can be neglected on orbital timescales. Note that k, which measures 
the strength of the tide, drops out. Now it is easy to show tha1@ (rr) = |(I — hh), which, in 
combination with (|22p . yields fl p -u+(ft p -h)(u • h) = 0. Solving for the equilibrium spin frequency 
of the planet, we have the correct statement of Cassini's first law: 

2n 

Q p = J^L (23) 

p cos (9 + sec 6» v ; 

This is the exact solution of the equations of motion when f2 p = 0. No equilibrium exists for 
9 > it/2; in particular, Cassini state 3 becomes unstable in the presence of the dissipative torque. 
Generally the synchronization time is comparable to the parallelization time (as we show next), so 
fl p — ► n as 6 — > 0, but if there is a mechanism for maintaining non-zero obliquity, like trapping 
in Cassini state 2, Q p < n. Equation (|23p can be understood as pseudo-synchronization when 
obliquity is held at a constant value; an obl ique and tidall y evolv ed planet can be expected to be 



a slow rotator. These points were made by iLevrard et al.l (120071 ). who gave an expression for 



with non-zero eccentricity (the dependence on 6 is separable from the dependence on e). 



Neglecting orbital precession, the dissipative torque from equations (|T5|) - (|T7|) will cause par- 
allelization; the component of the spin perpendicular to the orbit (£l p ±) decays approximately 
exponentially on a timescale t± = 2C p n(fih)~ 1 tp p . The deviation from synchronization {Q p h — n) 
also decays, but on a timescale t s = \t±. The twice longer timescale for parallelization is because 
the parallelization torque takes on smaller values in parts of the orbit, whereas the synchronization 
torque is constant. For precessing orbits that trap the spin in Cassini state 2, equation (|23p applies 
after a few parallelization times rj_, although it has no dependence on the tidal dissipation rate. 



4.2. Cassini states with dissipation 

The Cassini state formalism may be generalized to include ongoing tidal torques. For the 
generalized Cassini state we require that all the angular momenta vectors are stationary with 
respect to each other, even as the whole system precesses uniformly about the J axis. Cassini 



4 By symmetry, (ff) annihilates h and is a multiple of the identity I for any vector lying in the orbital plane. The 
factor (1/2) follows by considering the orbit average of (d ■ r) 2 for any constant vector d in the plane. 
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states are denned on timescales over which the orbital precession frequency g is nearly constant; 
on tidal dissipation timescales, g varies and the state evolves. 

In the presence of dissipation, fi p is shifted from the usual Cassini state; it does not lie in the 
plane defined by h and fi*. The reason is that dissipation introduces terms that attempt to align 
the planet's spin with its orbit, which would cause f2 p to slowly spiral in towards h. In particular, 
the second terms of equations (|16[) and (|17p cause a damping of the e- and q-components of the 
spin, according to equation (|15p . The system is able to avoid this drift by readjusting to a state in 
which the planetary spin executes pure precession about J. If the planet's spin and the orbit normal 
are stationary in a frame precessing with uniform angular velocity g = gj, then Ct p = g X fi p and 
h = g X h. Three constraints are needed on the three components of f2 p . 

To analyze the situation, let us align the reference direction q such that J lies in the q — h 
plane, making an angle I with the orbit normal h as in Figure Q] (i.e. q = Y an d e = X). 
Anticipating that Cassini state 2 will be shifted out of the q — h plane, denote the phase shift angle 
(f> s = 7r — cj), as it is the supplement of 4>. Let the angle between Ct p and h be denoted 9 and the 
angle between fi* and h be denoted 6*. Then, in the eqh coordinate system, we have 

g = ff(0, sin/, cos I), (24) 
ftp = Qp(sm <j) s sin#, cos 4> s sin 9, cos 9), (25) 
ftp = gflp (cos 4> s sin 9 cos I — cos 9 sin/, — sin</> s sin 9 cos /, sin (j) s sin 9 sin/). (26) 

Note that g < because the nodes regress. Equations (fl~5|) . (fTBp . and (fT7|) give the ratio between 
(f2 p ) e and {Ctp) q . Equating that to the ratio given by (|26p . we obtain the first constraint: 

— cos 4> s cos 9 + £ -1 sin cj) s — cos 4> s sin 9 cos / + cos 9 sin / 

smcj) s cos 9 + £ -1 cos (f) s sin <p s sin 9 cos / 

or equivalently, sin(c) + cj) s )/sm(5) = tan 9 cot/, (27) 

where 

M { Rn\^ 2 fl M 

^ = 2k p fl p tFp I — 1 = -Q p — — - and tan 5 = (£ cos 9)~ l . 

M p \ a J 3 n fi 

Note that £ has been written in terms of Q p rather than Q' p = 3Q P / '{Ak p ); nevertheless, it is likely 
to be quite large. 

The second constraint, similarly obtained from the ratio between (Ctp)^ and (Ct p ) q , gives an 
expression for the equilibrium spin rate that we have already derived in equation (|23|) . The final 
constraint can be derived by setting equal the precession rates of the planet's orbit and spin. 
The former (from Eq. [T4"|) is essentially due to the stellar oblateness since the rotational angular 
momentum of the planet is negligible {X p <C etc.): 

1 |h| 

9 = 



sin / I hi 



M+ V a ) n sin/ 
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The planet's spin precesses about the total angular momentum, maintaining with it an angle 
(3 = cos -1 (J • Cip), at the rate 



1 |»p| 
sin/3 |fl p | 

fj,h sin6 f r M / R p \ 5 



C p sin (3 \ L M p \ a J 1 + sec 



2 J\if 
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i 



2nt cos l 



1/2 



(29) 



Solutions of equations ([23]) . (|27|) . ([28]) . and ([29]) are generalized Cassini states. 



Given the complexity, solutions for 4> s and will be found by making a few approximations. 
To first order in <f) s and 5, which are small when dissipation is weak, (I27p becomes 

<p s w (Ccos6»)~ 1 (tan6»cot/- 1). (30) 

Setting equal (|28p and (|29p with the above approximation gives the following constraint for 0: 
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(31) 



In this final equation, the right hand side is determined by parameters that are approximately 
constant over a planetary spin precession period. On the left hand side, the quantities that are 
uncertain are the planet's apsidal-motion constant and its tidal friction timescale. If the former 
is too big, or the latter is too small, no real solution for e is possible, and there is no equilibrium 
corresponding to Cassini state 2. Tidal friction plays an important role here because with tp p — ► oo, 
there is always a solution for e. Indeed, equation © is recovered in this limilH. 



Starting from the values for Q p , 4> s , and e given by equations (|23|) . (|30l) and (|3Tj) . a Newton- 
Rap hson solver can efficiently find any solution to equations (|23l) . (f27l) . (f28l) . and (f29j) . Once e is 
found, (|23|) can be used to find the planetary spin rate. Our derivation gives an especially physical 
understanding for why the spin cannot be truly synchronous: there needs to be a small component 
of dp in the direction of the orbit normal so that the spin performs precession about the total 
angular momentum rather than the orbit normal. 



5 The spin precession constant a is related to the apsidal motion constant by 

fih M ( R P Y Q p 
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5. Evolution of Cassini states 

First we will show the integration of equations (|14p and ()15[) for an example that has an 
extremely tilted Cassini state 2. To find such a state we needed to take an extreme value for 
k p = 5 x 10 -5 , which corresponds to a polytrope n = 4.65, instead of a more reasonable value such 
as k p = 0.26 for n = 1. This choice lessens the planetary spin precession torque and dissipative 
torque relative to the rather weak orbital precession torque arising from the rotational bulge of the 
star. We also chose a quite large value of Q' p = 2.3 x 10 7 . The planet started with a semimajor 
axis of a = 0.05 AU. The stellar spin began with a period of 10 d and made an angle 45° with the 
orbit normal. The stellar apsidal motion constant was set to k+ = 0.0144, appropriate for an n = 3 
polytrope. The masses were M* = M & and M p = 1O _3 M0, radii i?* = 1.0R® and R p = O.IRq, 
and normalized moments of inertia C*/(M+R+) = 0.08 and C p /(M p Rp) = 0.25. 

See Figure 2] for the evolution into Cassini state 2. At first the planetary spin axis librates 
with large amplitude, tracing curves of nearly constant Hamiltonian (Eq. H|). After about 20 Myr, 
dissipation has caused the system to settle to a state of quasi-equilibrium (for comparison, t± = 4.3 
Myr). Panel (a) shows the evolution of the angles between the total angular momentum J and the 
orbit normal n, between J and the stellar spin axis fi*, and between J and the planetary spin axis 
fip. Panel (b) is a parametric plot (with time as the suppressed parameter) of the projection of 
the planetary spin vector f2 p onto the XY plane of Figure [TJ Panel (c) zooms in on the final state 
of panel (b), revealing that, in the presence of ongoing tidal dissipation, the equilibrium Cassini 
state 2 no longer lies on the y-axis; it is shifted towards positive X values by the tidal torque. 

The long-term consequence of such a state will be shown next. We started with the same initial 
conditions, except set Q' p = 4.6 x 10 6 so that the evolution would proceed faster. An integration for 
1.5 Gyr is shown in Figure Panels (a) and (b) are the same as those of Figured The timescale 
of damping into the Cassini state is a few Myr, smaller than the time axis resolution. Cassini 
state 2 now lies very far from the y-axis, and evolves even farther from it. The state "breaks" 
when Y = 0, then librates around, and damps into, Cassini state 1. The physical cause is that 
the tidal torque (the last terms on the right hand sides of Eqs. 1161 and 1171 which are proportional 
to a~®l 2 Q'~ l ) overcomes the orbital precession torque (the first terms on the right hand sides of 
Eqs. [TBI and W7\ which are proportional to a~ 3 ), as a decreases due to energy dissipation. Panel 
(c) shows the time dependence of the phase shift cf) s away from the X axis (c/. Fig. [TJ). Most of 
the time is spent in the evolving equilibrium of state 2, and the brief libration epochs that are 
prominent in panel (b) are not resolved in this plot. Panel (d) shows the semi-major axis evolution: 
dissipation in the planet converts orbital energy to heat energy. Panel (e) shows the evolution of 
the various angular momenta of the problem, projected along the total angular momentum axis. 
The angular momentum lost by the planetary orbit is absorbed by the stellar spin. However, in the 
absence of extremely strong dissipation in the star, the torque from the planet cannot cause the 
stellar spin period to change (Eq. [T3j) . Panel (f) shows a parametric plot of the stellar spin vector 
broken into two components: those parallel to and perpendicular to the planet's orbit. It is clear 
that the star's spin absorbs the angular momentum of the planetary orbit by reorienting, not by 
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changing magnitude. 

The non-zero value of <p s causes the stellar spin to come into alignment with the planet's orbit, 
without any dissipation in the star, by the following mechanism. The spin angular momentum 
of the planet is not coplanar with the other angular momenta, which causes the total angular 
momentum J to lie slightly outside of the plane defined by the dominant angular momenta L and 
S*. The evolution equations cause S* to precess about L (there is no spin-spin coupling), but 
this non-coplanarity directs a small component of that torque towards J. Similarly, the motion 
of the orbital angular momentum is not quite the pure precession of equation (|28p because the 
small coupling of the orbit to the rotational bulge of the planet (X p and y p of Eq. [T4"]) causes a 
secular motion of the orbit normal towards J. Therefore the stellar spin and planet orbit come 
into alignment with the total angular momentum, and thus with each other. Recall that the 
ingredients for this alignment are a close planet whose spin is trapped in Cassini state 2, whose 
tidal dissipation causes that state to shift in phase. In the example we have given, a significant spin- 
orbit misalignment persists because the planet's spin left Cassini state 2 before all of the available 
angular momentum from the star was tapped. 

In integrations like that of Figured] (in which <j) s = 0.1106), the first-order estimate of <p s 
(Eq. [30]) is quantitatively quite good, and its fidelity for larger values can be seen in Figure 02c). It 
is remarkable that the total amount of orbital energy dissipated is many orders of magnitude beyond 
that which is stored in the planet's spin; it is approximately the binding energy of the planet! The 
average dissipation rate indicated by the orbital migration of Figure 03(d) is E « 6.2 x 10 43 erg s^ 1 , 
enough to inflate the planet's radius (during the first 1.1 Gyr after Cassini state 2 is established). 
However, we achieved this result by adopting an apsidal moment constant almost four orders of 
magnitude smaller than Jupiter's, corresponding to a very strongly centrally concentrated body, in 
order to reduce the spin precession rate and dissipative tidal torque. Such a value for k p would be 
much more strongly discrepant with theory than the inflated radius itself. 

To seek the range of structural parameters that allow Cassini state 2 to be an equilibrium, we 
may use (j> s = 1 in equation (|30p to approximately indicate when Cassini state 2 will break due 
to dissipation in the planet. This first order estimate actually is serendipitously rather precise: 
recognizing from Figure [5] that Cassini state 2 actually breaks when <j) s = tt/2, we may substitute 
this value into equation ([2?]) . The result is that the right hand side of ([30]) equals 1. That is, the 
full phase shift equation gives <p s = tt/2 when the first order estimate gives <p s = 1, which can be 
seen in Figure 0(c) . 

What is the maximum energy dissipation rate a planet may experience while in Cassini state 2, 
driven by stellar bulge orbital precession? A closed equation for E max may be written by making 
a few approximations (all of which, we have verified, are well-justified). For realistic planetary 
structure, k p = 0.26, which requires 6 ~ tt/2. In the evaluations below, we will also use k± = 0.0144, 
R* = Rq, M* = M Q , CJ{M*Rl) = 0.08, SI* = 2^/(20d), 0* = 0.1, R p = Rj, M p = Mj, 
C p /(M p R 2 p ) = 0.25, and n = 2vr/(4d). Expanding g (Eq. [29]) to first order in tt/2 - 0, while 
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neglecting the dissipative term and setting sin/5 — > 1, then setting the resulting expression equal 
to equation (|28|) gives 

/ fc»itf aCpO^ sin 8* cos 0,, \ 1/2 
* ~ ff/2 - GM£ sin/ J ■ (32) 

Second, the solution will have an extreme value for (j) s , which we may identify by setting the first 
order estimate (Eq. l30l) to unity. This allows us to identify the Q p parameter for such a maximally 
dissipating planet: 

_3M, (vr/2-^3 
4 M tan 7 

For the nominal values above, this means Q p > 8 x 10 8 for Cassini state 2 to be an equilibrium, 
where equality holds for a state barely in equilibrium, with the planet dissipating at the maximal 
value. That the population of hot Jupiters is inferred to have Q p ~ 10 6 (£ )3.3p implies that a 
planet must be unusually non-dissipative to remain in state 2. Finally, the energy dissipation 
-Emax = G 2 M+M p Mh~ 3 h may be written 



^max = {V/m'^GMtRla^ 



n 



2 l/2 k 3/2 M 3/ 2 R 15/2 c 3/2^3 

/ sin 3 9+ cos 3 9+ \ 
I.V2 Ml o5/2 a 6 I gin /cos 2 / 



(34) 



which evaluates to E max 5 x 10 23 erg s^ 1 given the above values for the parameters. The 
maximum possible rate of heating due to orbital precession-stabilized obliquity tides is comparable 
to the luminosity of gravitational contraction. That is, if a certain transiting planet is in Cassini 
state 2, then tidal dissipation is not strong enough to inflate that planet's radius. Equation (|34p has 
a strong dependence on semimajor axis, but even a Jupiter mass planet at its Roche limit cannot 
dissipate a structurally significant power via Cassini state 2, even if it somehow managed to have 



Another way to assess stability is to follow lLevrard et al.l (|2007l ) in requiring that the dissipative 
torque does not destroy 9 stability (that a real solution exists for 6 in Eq. [3Tj) . For the above 
parameters, this gives E max ~6x 10 erg s , a factor of ~ 10 less restrictive than our criterion 
that the first order estimate of (f> s is less than 1. We interpret this to mean that Cassini equilibrium 
is generally lost in the (ft direction (as in Fig. E[b]), not the 9 direction. 



6. Discussion and Conclusions 

We have calculated the stability and evolution of Cassini state 2 with large obliquity 9 in the 
presence of tidal dissipation. We have verified that in principle, a spin state can couple the orbit 
to internal energy such that many times the energy of the spin is dissipated. However, it cannot 
dissipate the energy rapid enough for that heat to make a structural difference. Here we apply the 
concepts we have developed to the archetypal inflated planet, HD 209458b, and find that neither 
the stellar bulge nor a second planet can be responsible for a planet-inflating Cassini state. 
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According to iBodenheimer et al.l (|2003l ) , if the large radius of HD 209458b is due to internal 
heating, the inferred heating rate E is at least 4 x 10 26 erg s" 1 . Can this E be continually supplied 
by obliquity tides over the age of the planet, without the Cassini state breaking? 

In $5] we evaluated the hypothesis that the stellar bulge causes the orbital precession, finding 
that the tidal torque causes Cassini state 2 to fall out of equilibrium. When numbers appropriate 
to the HD 209458 system — including a plausible planetary apsidal-motion constant — are used in 
equation (|34p . the maximum possible heating rate is 3 x 10 24 erg s _1 , a factor of ~ 100 too small 
to explain the inflated radius. 

Independent of this dissipative torque constraint, we note that the usable angular momentum 
in the spin, as the star comes into alignment, is quite small given the observationally indicated 
stellar obliquity. The angular momentum that the stellar spin can absorb from the orbit before 
complete alignment, using equation (jHJ) with Li — > S* and i re i,i — ► is 

AL = L f -L = (L 2 + S* + 2LS+ cos0*) 1/2 — S* - L . (35) 

So for HD 209458, taking 0* = 0.1, we have AL « -2.1 x 10 46 g cm 2 s" 1 . The quoted E implies 
L ~ —1.9 x 10 31 g cm 2 s~ 2 . Then, if the angular momentum is lost at a steady rate, the system 
will become coplanar and the high obliquity state will end in only ~ 35 Myr, which is < 1% of the 



system's age (5 Gyr: ICody <fe Sasselovll2002l ). Under this hypothesis we are seeing the system at a 



very special time, just before it becomes coplanar. 

Alternatively, a second planet, rather than stellar obliquity, may be responsible for the orbital 
precession that has HD 209458b trapped in Cassini state 2 (see §3.2|) . Here we compile the con- 
straints and find that this hypothesis is also ruled out (Fig. [6]). Briefly, a putative second planet 
must be massive and distant enough to absorb the angular momentum of the dissipating planet, 
yet close enough so that its torque can compete with the dissipative torque, and these requirements 
would violate the radial velocity constraints. Each of these constraints is explained next. 

How much angular momentum must be transfered to the second planet? If we happen to be 
observing the system near the end of Cassini state 2, such that it is nearly coplanar now, the final 
angular momentum in the orbit of HD 209458b is Lf = L now = 1.26 x 10 49 g cm 2 s" 1 . We may 
combine the current inferred rate of orbital evolution with the expected scaling of L with L to 
find its primordial value ( §3. 2j) . Assuming Cassini state 2 was established at the inception of the 
system, the planet's primordial angular momentum was Li ~ 1.41 x 10 49 g cm 2 s _1 . Thus the second 
planet must have enough non-aligned angular momentum such that ALi = 1.5 x 10 g cm 2 s" 1 
is absorbed before coplanarity is established. If we had assumed that the Cassini state was only 
midway through its evolution, this requirement would be even more strict. 

Using these values, we have that the primordial value for i re i j must have been greater than 
i cr it = 27°. This value is considerably larger than relative inclinations of planets in our solar system, 
but it may be envisioned that planet-planet scattering, followed by planet-star tidal dissipation, 
was the cause of HD 209458b's current small semi-major axis, after which i m \ may plausibly be 



-23- 



that large (jMarzari Sz Weidenschilling||2002l ). On the other hand, for very large initial inclinations, 
eccentricity oscillations by the Kozai mechanism plus tidal friction will reduce i re \ to ~ 40° before 
obliquity tides start to dominate dissipation, so it is reasonable to consider only initial values of 
i re \ t i less than 40°. 

For each assumed primordial relative inclination, the second planet must have at least the 
orbital angular momentum L2 of equation ([9|) so that the two orbits are not yet coplanar. This 
lower limit gives a lower limit to the mass of the putative planet as a function of 02(1 — e?,), which 
we plot in Fig ure [6] for two values of i re i i> 

At large orbital separation, the torque of the second planet will not be able to compete with 
the dissipative torques, and by analogy with the stellar oblateness analysis, Cassini state 2 will 
not be an equilibrium. Here equation (|30|) still applies, but 8 must be determined by setting the 
planet's spin precession rate g of equation (j29j) equal to 

_ 3 G 1 / 2 M 2 a 3 / 2 coBi re i, now 
92 4 MV2 a 3 (l-e 2 ) 3 / 2 ' {6b) 

which is the orbital precession rate due to the second planet, whose mass is M2, semimajor axis 
02, eccentricity e2, and current relative inclination i re i,now Here 52 is calculated by secular theory, 
and any enhancement to this rate by the effects of mean motion resonance will only be temporary, 
since a is evolving due to tidal dissipation. To find the critical parameters of a second planet 
which barely compete with the dissipation torque, we may set (fi s = 1 in equation (|30p . and similar 
approximations as those that led to equations (|32p - (|34p give the constraint 

M 2 > 2x4y^ 2/3 M 1 J 3 k 1 p / 3 R 5 / i cot 2 / 3 / ^ 



a^l — e 2 ,) 3 / 2 '3 G C p a coszre^. 

If the exterior planet's orbit dominates the angular momentum, i ro i,now ~ I] the minimum value of 
cot 2 / 3 I /cos I occurs at / = tan^ 1 ^2 « 54.7° and has the value 3 1/2 /2 1 / 3 « 1.375. This angle is 
bigger than the critical angle for Kozai eccentricity oscillations, so taking 1.375 for the multiplicative 
factor depending on the inclinations is quite conservative. Then, evaluating the other terms with 
the values cited for HD 209458b gives il^a^l - e|)" 3/2 > 3.74Af/AlT 3 . This constraint is 
labeled "dissipative torque" on Figure El It scales with e2 differently than the angular momentum 
constraint, but the limit shown is true for all eccentricities. 

The parameter space of possible second planets is also constrained by the lack of a secondary 
signal in the radial velocity data, giving a mass upper limit plotted in Figure [H No part of parameter 
space remains, so the second planet mechanism for causing the inflation of HD 209458b via Cassini 
state 2 is ruled out. 

We have ruled out both the stellar rotational bulge and a second planet as drivers of a high 
obliquity Cassini state for HD 209458b, if it is presumed that obliquity tides are inflating its radius. 
We have focused on this planet as the observational material far surpasses that of other planets. 
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However, we expect that the hypothesis also fails for the other anomalously large planets HAT- 
P-lb and WASP- lb. Likewise, planetary systems discovered in the future will need very peculiar 
parameters if obliquity tides are to be an important heating source. 
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Fig. 1. — The XYZ coordinate system precesses with angular frequency g about the direction of 
the total angular momentum J. The orbital plane, with normal h along the Z-axis, is inclined to 
the invariable plane by angle I with the ascending node in the X direction (J lies in the Y Z plane) . 
The xyz axes are the principal axes of the planet, with the unit spin vector S along the axis of 
symmetry. (j>, if), and the obliquity 9 are Euler angles. The phase angle <p s = tt — (j). Nondissipative 
Cassini states have S in the YZ plane, with <f) s = tt (0) for state 1 (2); but dissipation causes 
< (f> s < tt/2 in state 2. 
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Fig. 2. — Obliquities of the four Cassini states as functions of —g/a. Solid numbered curves are 
the solutions of equation Q for / = 0.1 (left panel) and I = 1.0 (right panel). Cassini states lie in 
the Y Z plane of Figure [TJ by convention 9 is positive when <p s = and negative when <j) s = ir. The 
separatrix intersects the YZ plane at state 4, and at the angles indicated by the dashed curves. 
States 1 and 4 merge and vanish at the critical value (— g/a) CT it = (sin 2 / 3 / + cos 2 / 3 J)~ 3 / 2 (dotted 
vertical line). Three values of 9 are indicated by dotted horizontal lines. The largest absolute 
obliquity attained by state 1 is 9 cr \t — tan _1 (— tan 1 / 3 /). States 2 and 3 approach 9 = 1 and 
9 = I — ir, respectively, when \g/a\ 3> 1- 
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Fig. 3. — Measured stellar rotational velocity (usini)* relative to synchronous rotation v sy nch = 
nR-k (squares). Stars bearing transiting planets are listed top to bottom in order of decreasing 
v sinz/'u sync } 1 . Error bars indicate observational uncertainty in (usini)* and R*. We compute 

^brake from a Skumanicll law (see text) and plot fbrake/^synch 

(crosses). Scatter in f brake is ~ 
10%. Gray boxes give a range for Wbrake/^synch corresponding to the 1 a uncertainty in stellar age 
(systematic uncertainties are not accounted for). Probably all stars are rotating subsynchronously, 
since v sin i k, v for a transiting planet whose orbit is nearly coplanar with the equator of its host 
star. Only HD 209458 is rotating ( marginally) more rapidl y than would be predicted by magnetic 
braking alone. Data are taken from lBurrows et al.l (|2006al ) and references the rein, except for th e 
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Fig. 4. — Settling into Cassini state 2. (a) The angles between J and the other vectors as a function 
of time, (b) The motion of f2 p on the unit sphere, in the coordinate system of Figure [TJ (c) After 
the equilibrium state is reached, the system slowly evolves (see Fig. [5]). Notice that the state does 
not have equilibrium value (f> = tv, as expected from the analysis with no dissipation. In this case 
there is a phase shift away from the Y axis of <p s = 0.1105. We have used an exceptionally low 
value of the apsidal motion constant k p = 5 x 10 -5 for this integration (see text). 
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Fig. 5. — The long term evolution of Cassini state 2. (a) The angles between J and the other 
vectors as a function of time, (b) The motion of Cl p on the unit sphere, in the coordinate system of 
Figure [TJ This longer integration shows that after damping into Cassini state 2, as in Figure IU(b), 
the planetary spin vector evolves to negative Y, after which equilibrium is lost it librates around, 
and damps into, Cassini state 1. (c) The evolution of the phase shift angle (f) s = ir — (f>. Most of 
the time in the parametric plot of panel (b) is spent in the evolving state 2; libration epochs are 
too short to be resolved on this plot. The dashed line is the first-order estimate of equation (I30p . 
(d) The obliquity tides cause the semi-major axis to shrink. After damping to Cassini state 1, no 
appreciable energy is lost, (e) The changes in the projections of the three angular momenta along 
the J axis, as a function of time. The angular momentum lost by the planet as it migrates in is 
absorbed by the stellar spin, (f) The stellar spin absorbs the angular momentum of the planetary 
orbit by reorienting toward alignment. The dashed line is the set of points with equal magnitude of 
0*, showing that the only change to the stellar spin is its direction. We have used an exceptionally 
low value of the apsidal motion constant k p = 5 x 10 -5 for this integration (see text). 
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Fig. 6. — The limits of mass M2 and orbital elements 02(1 — e|) of a second planet for the 
hypothesis that it is responsible for the inflation of HD 209458b via a Cassini resonance. Here i ie n 
is the initial relative inclination between the two planets, and the bounds so labeled reflect the 
tendency of the second planet's orbit to align as it absorbs the angular momentum shed by the 
orbit of HD 209458b (L2 must be greater than that of eq. [9]). Values of i re i,i > 40° are unlikely 
as the second planet would force Kozai eccentricity oscillations in the transiting planet. The curve 
labeled "dissipative torque" is the lower limit of mass needed to overcome the tidal dissipation 
that tends to break Cassini state 2 (eq. [37]). The radial velocity dataset also constrains possible 
second planets; absence of a secondary signal of amplitude K = 5 m s _1 puts an upper limit on 
mass at each orbital distance. The arrows point to the portion of parameter space allowed by each 
particular constraint. Apparently all of parameter space is ruled out: tides due to a high obliquity 
state, stabilized by the orbital precession induced by a second planet, cannot significantly inflate 
HD 209458b. See ^6] for derivation and details. 



